This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.
Issue Howard Bayes
Time series 1999 - present 1977 - present
Bias in SWHS Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests. Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data. Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas. Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with. By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.
**Figure 1.**- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.


**Figure 2.**- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).


Results

**Figure X.**- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

**Figure 3.**- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 3.**- Total rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

**Figure 4.**- Black rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


And black rockfish releases…

**Figure 5.**- Black rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 6.**- Yellow rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 7.**- Yellow rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 8.**- DSR rockfish (excluding yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 9.**- DSR rockfish releases (including yelloweye) 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 11.**- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 12.**- Slope rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

**Figure 13.**- *Other* rockfish (non-pelagic and non-yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 14.**- *Other* rockfish (non-pelagic and non-yelloweye) releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



Total Biomass Removal Estimates

**Figure 15.**- Black rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.



**Figure 16.**- Yellow rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

**Figure 17.**- Pelagic rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 18.**- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 19.**- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


Model fit

Logbook residuals

**Figure 20.**- Residuals from logbook harvests.

Figure 20.- Residuals from logbook harvests.


SWHS residuals

**Figure 21.**- Residuals from SWHS harvests.

Figure 21.- Residuals from SWHS harvests.



**Figure 22.**- Residual of SWHS releases.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

**Figure 23.**- Mean percent of harvest by charter anglers.

Figure 23.- Mean percent of harvest by charter anglers.


When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

**Figure 24.**- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.


**Figure 25.**- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 26.**- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 27.**- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.


**Figure 28.**- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.


**Figure 29.**- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

**Figure 30.**- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.


Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

**Figure 31.**- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

**Figure 32.**- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

**Figure 33.**- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

**Figure 34.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

**Figure 35.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

**Figure 36.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.



P(slope|non-pelagic & non-yellowye) For release estimates

## NULL



Weight Fits

**Figure 37.**- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 38.**- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 39.**- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 40.**- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 41.**- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.4692982

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.06 while 0.94 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable

area_n

year

user

mean

median

lower

upper

Rhat

sd

skew

area

p_yellow

11

27

2

0.435

0.423

0.175

0.748

1.013

0.150

0.080

CSEO

p_pelagic

13

1

1

0.268

0.243

0.065

0.598

1.199

0.143

0.175

NSEI

p_pelagic

13

1

2

0.108

0.088

0.022

0.303

1.176

0.075

0.265

NSEI

p_yellow

13

1

1

0.590

0.543

0.158

0.993

1.021

0.280

0.168

NSEI

p_yellow

13

1

2

0.590

0.539

0.170

0.993

1.019

0.278

0.183

NSEI

p_pelagic

13

10

1

0.291

0.270

0.085

0.592

1.144

0.136

0.154

NSEI

p_pelagic

13

10

2

0.119

0.102

0.028

0.305

1.163

0.075

0.228

NSEI

p_yellow

13

10

1

0.592

0.546

0.150

0.992

1.021

0.279

0.165

NSEI

p_yellow

13

10

2

0.591

0.554

0.146

0.993

1.019

0.282

0.131

NSEI

p_pelagic

13

11

1

0.297

0.281

0.088

0.588

1.143

0.132

0.121

NSEI

p_pelagic

13

11

2

0.121

0.103

0.029

0.311

1.127

0.073

0.246

NSEI

p_yellow

13

11

1

0.594

0.549

0.154

0.993

1.019

0.280

0.161

NSEI

p_yellow

13

11

2

0.593

0.551

0.158

0.993

1.020

0.279

0.150

NSEI

p_pelagic

13

12

1

0.300

0.282

0.091

0.600

1.146

0.136

0.133

NSEI

p_pelagic

13

12

2

0.123

0.107

0.029

0.304

1.150

0.072

0.222

NSEI

p_yellow

13

12

1

0.590

0.541

0.165

0.993

1.019

0.280

0.175

NSEI

p_yellow

13

12

2

0.590

0.540

0.156

0.993

1.020

0.281

0.178

NSEI

p_pelagic

13

13

1

0.306

0.290

0.101

0.603

1.117

0.133

0.121

NSEI

p_pelagic

13

13

2

0.126

0.110

0.032

0.313

1.132

0.074

0.215

NSEI

p_yellow

13

13

1

0.591

0.547

0.153

0.992

1.021

0.280

0.157

NSEI

p_yellow

13

13

2

0.590

0.538

0.148

0.993

1.023

0.281

0.185

NSEI

p_pelagic

13

14

1

0.291

0.273

0.098

0.562

1.121

0.123

0.147

NSEI

p_pelagic

13

14

2

0.131

0.114

0.033

0.319

1.125

0.074

0.229

NSEI

p_yellow

13

14

1

0.589

0.545

0.163

0.993

1.022

0.279

0.158

NSEI

p_yellow

13

14

2

0.592

0.549

0.153

0.994

1.024

0.280

0.153

NSEI

p_pelagic

13

15

1

0.270

0.255

0.092

0.541

1.111

0.118

0.127

NSEI

p_pelagic

13

15

2

0.120

0.106

0.035

0.283

1.125

0.065

0.216

NSEI

p_yellow

13

15

1

0.590

0.547

0.160

0.992

1.022

0.280

0.154

NSEI

p_yellow

13

15

2

0.588

0.544

0.158

0.993

1.019

0.280

0.157

NSEI

p_pelagic

13

16

1

0.309

0.298

0.116

0.571

1.081

0.121

0.091

NSEI

p_pelagic

13

16

2

0.138

0.122

0.039

0.319

1.106

0.074

0.216

NSEI

p_yellow

13

16

1

0.588

0.546

0.159

0.993

1.020

0.279

0.151

NSEI

p_yellow

13

16

2

0.590

0.536

0.156

0.993

1.018

0.280

0.193

NSEI

p_pelagic

13

17

1

0.323

0.311

0.129

0.593

1.064

0.121

0.099

NSEI

p_pelagic

13

17

2

0.143

0.129

0.043

0.336

1.105

0.075

0.186

NSEI

p_yellow

13

17

1

0.592

0.553

0.157

0.992

1.024

0.281

0.139

NSEI

p_yellow

13

17

2

0.589

0.541

0.156

0.994

1.021

0.281

0.171

NSEI

p_pelagic

13

18

1

0.331

0.318

0.132

0.595

1.080

0.119

0.109

NSEI

p_pelagic

13

18

2

0.147

0.133

0.044

0.324

1.100

0.075

0.187

NSEI

p_yellow

13

18

1

0.593

0.556

0.163

0.993

1.021

0.280

0.132

NSEI

p_yellow

13

18

2

0.591

0.545

0.148

0.992

1.019

0.280

0.165

NSEI

p_pelagic

13

19

1

0.330

0.319

0.131

0.592

1.066

0.120

0.092

NSEI

p_pelagic

13

19

2

0.152

0.137

0.050

0.327

1.051

0.073

0.205

NSEI

p_yellow

13

19

1

0.588

0.543

0.157

0.993

1.021

0.281

0.160

NSEI

p_yellow

13

19

2

0.589

0.538

0.156

0.993

1.023

0.281

0.182

NSEI

p_pelagic

13

2

1

0.266

0.244

0.067

0.575

1.197

0.140

0.158

NSEI

p_pelagic

13

2

2

0.109

0.090

0.021

0.290

1.193

0.075

0.254

NSEI

p_yellow

13

2

1

0.590

0.549

0.156

0.993

1.017

0.280

0.146

NSEI

p_yellow

13

2

2

0.592

0.550

0.158

0.993

1.021

0.280

0.150

NSEI

p_pelagic

13

20

1

0.351

0.340

0.145

0.613

1.059

0.120

0.091

NSEI

p_pelagic

13

20

2

0.158

0.145

0.055

0.330

1.043

0.072

0.181

NSEI

p_yellow

13

20

1

0.589

0.537

0.157

0.992

1.018

0.281

0.185

NSEI

p_yellow

13

20

2

0.591

0.542

0.162

0.993

1.024

0.279

0.176

NSEI

p_pelagic

13

21

1

0.371

0.361

0.169

0.624

1.044

0.120

0.084

NSEI

p_pelagic

13

21

2

0.170

0.158

0.060

0.348

1.028

0.076

0.159

NSEI

p_yellow

13

21

1

0.587

0.537

0.164

0.993

1.028

0.279

0.179

NSEI

p_yellow

13

21

2

0.592

0.552

0.160

0.993

1.023

0.281

0.142

NSEI

p_pelagic

13

22

2

0.170

0.157

0.061

0.349

1.026

0.075

0.173

NSEI

p_yellow

13

22

1

0.586

0.539

0.151

0.993

1.022

0.282

0.166

NSEI

p_yellow

13

22

2

0.589

0.550

0.164

0.993

1.019

0.279

0.140

NSEI

p_yellow

13

23

1

0.562

0.514

0.125

0.991

1.022

0.290

0.166

NSEI

p_yellow

13

23

2

0.628

0.609

0.184

0.995

1.017

0.267

0.071

NSEI

p_yellow

13

24

1

0.565

0.514

0.131

0.990

1.026

0.287

0.178

NSEI

p_yellow

13

24

2

0.627

0.606

0.184

0.995

1.024

0.264

0.079

NSEI

p_yellow

13

25

1

0.567

0.523

0.131

0.991

1.024

0.284

0.155

NSEI

p_yellow

13

25

2

0.626

0.597

0.191

0.995

1.023

0.263

0.110

NSEI

p_yellow

13

26

1

0.554

0.511

0.129

0.991

1.027

0.287

0.150

NSEI

p_yellow

13

26

2

0.624

0.599

0.189

0.994

1.018

0.263

0.095

NSEI

p_yellow

13

27

1

0.538

0.482

0.117

0.990

1.031

0.286

0.196

NSEI

p_yellow

13

27

2

0.614

0.581

0.183

0.994

1.025

0.261

0.126

NSEI

p_yellow

13

28

1

0.545

0.501

0.128

0.989

1.032

0.272

0.162

NSEI

p_yellow

13

28

2

0.610

0.584

0.185

0.993

1.028

0.255

0.102

NSEI

p_yellow

13

29

1

0.499

0.444

0.126

0.983

1.048

0.256

0.215

NSEI

p_yellow

13

29

2

0.581

0.553

0.182

0.990

1.035

0.244

0.115

NSEI

p_pelagic

13

3

1

0.271

0.250

0.074

0.582

1.137

0.136

0.154

NSEI

p_pelagic

13

3

2

0.110

0.091

0.022

0.297

1.181

0.073

0.261

NSEI

p_yellow

13

3

1

0.591

0.550

0.150

0.993

1.019

0.281

0.146

NSEI

p_yellow

13

3

2

0.595

0.555

0.153

0.993

1.019

0.279

0.143

NSEI

p_pelagic

13

4

1

0.271

0.249

0.069

0.588

1.188

0.139

0.159

NSEI

p_pelagic

13

4

2

0.111

0.091

0.023

0.298

1.159

0.074

0.269

NSEI

p_yellow

13

4

1

0.588

0.545

0.156

0.993

1.021

0.282

0.153

NSEI

p_yellow

13

4

2

0.590

0.546

0.149

0.993

1.021

0.283

0.156

NSEI

p_pelagic

13

5

1

0.277

0.253

0.072

0.594

1.186

0.139

0.173

NSEI

p_pelagic

13

5

2

0.111

0.094

0.023

0.302

1.166

0.072

0.235

NSEI

p_yellow

13

5

1

0.592

0.559

0.159

0.992

1.023

0.279

0.118

NSEI

p_yellow

13

5

2

0.594

0.547

0.161

0.993

1.022

0.279

0.168

NSEI

p_pelagic

13

6

1

0.276

0.254

0.074

0.589

1.165

0.138

0.159

NSEI

p_pelagic

13

6

2

0.112

0.093

0.023

0.293

1.172

0.073

0.261

NSEI

p_yellow

13

6

1

0.591

0.549

0.150

0.993

1.019

0.280

0.150

NSEI

p_yellow

13

6

2

0.592

0.547

0.158

0.993

1.021

0.280

0.161

NSEI

p_pelagic

13

7

1

0.280

0.261

0.075

0.582

1.205

0.138

0.138

NSEI

p_pelagic

13

7

2

0.115

0.098

0.025

0.301

1.148

0.073

0.232

NSEI

p_yellow

13

7

1

0.590

0.547

0.149

0.993

1.022

0.280

0.153

NSEI

p_yellow

13

7

2

0.596

0.552

0.163

0.993

1.021

0.278

0.158

NSEI

p_pelagic

13

8

1

0.284

0.267

0.081

0.584

1.156

0.136

0.125

NSEI

p_pelagic

13

8

2

0.116

0.097

0.026

0.307

1.172

0.075

0.254

NSEI

p_yellow

13

8

1

0.592

0.553

0.158

0.992

1.019

0.280

0.139

NSEI

p_yellow

13

8

2

0.595

0.551

0.157

0.993

1.020

0.279

0.157

NSEI

p_pelagic

13

9

1

0.290

0.271

0.083

0.592

1.159

0.138

0.137

NSEI

p_pelagic

13

9

2

0.119

0.101

0.026

0.309

1.180

0.074

0.242

NSEI

p_yellow

13

9

1

0.589

0.546

0.153

0.993

1.023

0.281

0.153

NSEI

p_yellow

13

9

2

0.590

0.553

0.158

0.993

1.025

0.281

0.132

NSEI

p_yellow

14

7

1

0.575

0.563

0.259

0.966

1.010

0.181

0.066

NSEO

p_pelagic

16

1

1

0.390

0.386

0.152

0.669

1.015

0.136

0.029

SSEO

p_pelagic

16

1

2

0.370

0.361

0.144

0.639

1.013

0.133

0.068

SSEO

p_slope

16

10

2

0.079

0.058

0.010

0.264

1.013

0.070

0.301

SSEO

p_pelagic

16

14

2

0.374

0.365

0.150

0.638

1.012

0.128

0.070

SSEO

p_pelagic

16

3

1

0.387

0.382

0.146

0.663

1.014

0.136

0.037

SSEO

p_pelagic

16

3

2

0.367

0.354

0.134

0.638

1.010

0.133

0.098

SSEO

p_pelagic

16

4

2

0.368

0.356

0.141

0.650

1.011

0.134

0.090

SSEO

p_pelagic

16

5

1

0.392

0.389

0.150

0.664

1.012

0.135

0.022

SSEO

p_pelagic

16

5

2

0.370

0.363

0.142

0.639

1.011

0.132

0.053

SSEO

p_pelagic

16

6

1

0.390

0.386

0.154

0.657

1.011

0.133

0.030

SSEO

p_pelagic

16

9

1

0.394

0.386

0.158

0.665

1.011

0.134

0.060

SSEO

p_pelagic

3

1

1

0.541

0.546

0.267

0.789

1.013

0.137

-0.036

PWSI

p_pelagic

3

1

2

0.184

0.170

0.059

0.373

1.011

0.084

0.166

PWSI

p_pelagic

3

12

1

0.543

0.548

0.280

0.782

1.010

0.132

-0.038

PWSI

p_pelagic

3

16

2

0.185

0.173

0.062

0.384

1.012

0.083

0.144

PWSI

p_pelagic

3

3

2

0.186

0.172

0.062

0.388

1.013

0.086

0.163

PWSI

p_pelagic

3

6

2

0.183

0.171

0.060

0.370

1.018

0.081

0.147

PWSI

p_black

8

40

1

0.486

0.417

0.343

0.964

1.030

0.189

0.366

afognak

Table Xb. True parameters with an Rhat in excess of 1.1

variable

area_n

year/species

mean

median

lower

upper

Rhat

sd

skew

year

area

beta0_pH

1

1

0.584

0.591

0.212

0.910

1.117

0.178

-0.039

CI

beta1_pH

1

1

2.456

2.384

1.755

3.640

1.196

0.463

0.155

CI

beta2_pH

10

3

0.289

0.225

0.097

0.795

1.142

0.301

0.213

northeast

beta1_pH

11

1

3.363

3.275

2.619

4.768

1.120

0.522

0.169

CSEO

beta0_pelagic

13

-1.282

-1.266

-2.553

-0.140

1.318

0.678

-0.024

NSEI

beta2_pelagic

13

0.159

0.139

0.075

0.376

1.129

0.078

0.255

NSEI

beta3_pelagic

13

31.044

30.544

23.503

40.660

1.354

4.891

0.102

NSEI

beta2_yellow

13

-4.831

-3.201

-16.709

-0.318

1.172

4.687

-0.348

NSEI

beta2_yellow

14

-6.235

-4.935

-19.024

-0.354

1.109

5.092

-0.255

NSEO

beta0_pH

8

1

-0.654

-0.633

-1.430

-0.076

1.128

0.341

-0.062

afognak

beta0_H

8

5.922

6.467

-18.495

29.430

1.111

34.588

-0.016

afognak

beta2_pH

9

1

1.110

0.477

0.138

7.840

1.229

1.790

0.354

eastside

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.01

527.00

210.00

1,034.00

356.00

Prop. of total

0.72

0.28

0.74

0.26

Prop. median less than 1

0.25

0.21

0.30

0.43

Prop. median less than 10

0.35

0.56

0.50

0.74

Prop. median less than 50

0.46

0.80

0.62

0.81

Median of the median values

101.83

8.48

10.57

1.91

Mean of the median values

606.48

159.90

298.96

148.08

Mean skew of posteriors

0.23

0.29

0.27

0.25

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.1

39.00

15.00

57.00

40.00

Prop. of total

0.72

0.28

0.59

0.41

Prop. median less than 1

0.08

0.40

0.75

0.55

Prop. median less than 10

0.08

0.80

0.82

0.80

Prop. median less than 50

0.08

1.00

0.82

0.90

Median of the median values

338.91

2.11

0.04

0.55

Mean of the median values

629.04

8.02

46.97

8.96

Mean skew of posteriors

0.25

0.20

0.18

0.18

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable

Freq

beta0_H

1

beta0_pelagic

1

beta0_pH

2

beta1_pH

2

beta2_pelagic

1

beta2_pH

2

beta2_yellow

1

beta3_pelagic

1

Bp_ay

10

Bp_ayg

3

Bp_ayu

8

eps_pel

2

eps_pH

1

H_ayg

1

Hb_ay

10

Hb_ayg

6

Hb_ayu

7

Ho_ay

2

Ho_ayg

1

Ho_ayu

4

Hp_ay

15

Hp_ayg

9

Hp_ayu

11

Hs_ay

1

Hs_ayg

1

Hs_ayu

1

Hy_ay

4

Hy_ayg

1

Hy_ayu

6

p_pelagic

30

p_slope

1

R_ayg

4

R_ayu

2

Rb_ay

1

Rb_ayg

3

Rb_ayu

1

Rd_ayg

2

Rd_ayu

3

re_pelagic

8

re_pH

5

Ro_ay

13

Ro_ayg

7

Ro_ayu

30

Rp_ay

1

Rp_ayg

4

Rp_ayu

2

Rs_ayg

4

Rs_ayu

1

Ry_ay

12

Ry_ayg

4

Ry_ayu

26

Parameter estimates:

Summary Table of Parameter Estimates
Parameter mean sd Lower_CI Median Upper_CI
mu_bc_H[1] -0.019 0.069 -0.152 -0.022 0.122
mu_bc_H[2] 0.076 0.054 -0.024 0.073 0.192
mu_bc_H[3] -0.272 0.068 -0.399 -0.273 -0.139
mu_bc_H[4] -0.918 0.179 -1.278 -0.915 -0.570
mu_bc_H[5] 1.321 1.093 0.110 1.062 4.052
mu_bc_H[6] -2.425 0.314 -3.031 -2.431 -1.788
mu_bc_H[7] -0.455 0.103 -0.664 -0.452 -0.252
mu_bc_H[8] 0.245 0.350 -0.343 0.213 1.000
mu_bc_H[9] -0.243 0.126 -0.485 -0.243 0.014
mu_bc_H[10] 0.061 0.065 -0.058 0.058 0.194
mu_bc_H[11] -0.124 0.038 -0.198 -0.124 -0.049
mu_bc_H[12] -0.249 0.108 -0.479 -0.244 -0.048
mu_bc_H[13] -0.156 0.073 -0.299 -0.157 -0.015
mu_bc_H[14] -0.307 0.090 -0.489 -0.307 -0.140
mu_bc_H[15] -0.330 0.046 -0.418 -0.329 -0.236
mu_bc_H[16] -0.158 0.369 -0.783 -0.188 0.660
mu_bc_R[1] 1.040 0.120 0.809 1.040 1.269
mu_bc_R[2] 1.015 0.077 0.856 1.016 1.157
mu_bc_R[3] 0.221 0.091 0.041 0.225 0.393
mu_bc_R[4] 0.661 0.203 0.234 0.671 1.059
mu_bc_R[5] 1.674 0.366 0.946 1.680 2.389
mu_bc_R[6] -1.447 0.338 -2.099 -1.453 -0.772
mu_bc_R[7] 0.516 0.166 0.177 0.519 0.836
mu_bc_R[8] 0.644 0.150 0.362 0.639 0.939
mu_bc_R[9] 0.389 0.166 0.048 0.394 0.705
mu_bc_R[10] 0.897 0.101 0.701 0.895 1.089
mu_bc_R[11] 1.096 0.074 0.954 1.096 1.241
mu_bc_R[12] 0.874 0.161 0.564 0.871 1.187
mu_bc_R[13] 1.026 0.091 0.846 1.029 1.202
mu_bc_R[14] 0.852 0.121 0.607 0.856 1.080
mu_bc_R[15] 0.779 0.098 0.590 0.778 0.974
mu_bc_R[16] 1.262 0.111 1.043 1.261 1.477
tau_pH[1] 2.459 0.183 2.112 2.451 2.819
tau_pH[2] 1.564 0.148 1.295 1.556 1.873
tau_pH[3] 2.276 0.304 1.729 2.260 2.912
tau_pH[4] 1.114 0.175 0.810 1.100 1.496
tau_pH[5] 0.790 0.131 0.566 0.782 1.083
beta0_pH[1,1] 0.584 0.178 0.212 0.591 0.910
beta0_pH[2,1] 1.435 0.165 1.090 1.440 1.768
beta0_pH[3,1] 1.382 0.254 0.754 1.415 1.776
beta0_pH[4,1] 1.532 0.263 0.950 1.558 1.962
beta0_pH[5,1] -0.816 0.389 -1.756 -0.777 -0.177
beta0_pH[6,1] -0.560 0.453 -1.703 -0.454 0.042
beta0_pH[7,1] 0.752 0.010 0.732 0.753 0.772
beta0_pH[8,1] -0.654 0.341 -1.430 -0.633 -0.076
beta0_pH[9,1] -0.331 0.293 -0.984 -0.303 0.165
beta0_pH[10,1] 0.262 0.202 -0.172 0.276 0.622
beta0_pH[11,1] 0.043 0.322 -0.726 0.091 0.557
beta0_pH[12,1] 0.637 0.185 0.260 0.640 0.981
beta0_pH[13,1] -0.640 0.149 -0.957 -0.627 -0.368
beta0_pH[14,1] 0.006 0.172 -0.351 0.012 0.328
beta0_pH[15,1] 0.028 0.205 -0.442 0.049 0.371
beta0_pH[16,1] -0.860 0.223 -1.333 -0.846 -0.456
beta0_pH[1,2] 2.566 0.086 2.405 2.564 2.739
beta0_pH[2,2] 3.078 0.022 3.037 3.078 3.122
beta0_pH[3,2] 3.627 0.041 3.547 3.626 3.708
beta0_pH[4,2] 2.936 0.019 2.899 2.936 2.973
beta0_pH[5,2] 4.525 1.429 2.716 4.259 8.039
beta0_pH[6,2] 3.330 0.158 3.025 3.323 3.654
beta0_pH[7,2] 2.025 0.049 1.932 2.025 2.125
beta0_pH[8,2] 2.822 0.066 2.695 2.820 2.951
beta0_pH[9,2] 3.963 0.206 3.581 3.959 4.383
beta0_pH[10,2] 4.001 0.129 3.756 3.998 4.259
beta0_pH[11,2] 1.888 0.009 1.871 1.888 1.906
beta0_pH[12,2] 1.799 0.069 1.671 1.799 1.939
beta0_pH[13,2] 2.497 0.026 2.446 2.496 2.550
beta0_pH[14,2] 1.806 0.024 1.759 1.806 1.853
beta0_pH[15,2] 2.550 0.017 2.516 2.550 2.584
beta0_pH[16,2] 2.776 0.025 2.728 2.776 2.825
beta0_pH[1,3] 0.676 0.591 -0.897 0.830 1.373
beta0_pH[2,3] 2.126 0.018 2.092 2.127 2.160
beta0_pH[3,3] 2.689 0.024 2.641 2.689 2.735
beta0_pH[4,3] 2.981 0.035 2.916 2.980 3.049
beta0_pH[5,3] 3.828 2.670 0.988 3.100 11.191
beta0_pH[6,3] 1.253 0.326 0.683 1.229 1.955
beta0_pH[7,3] 0.855 0.050 0.758 0.854 0.954
beta0_pH[8,3] 0.284 0.119 0.032 0.290 0.502
beta0_pH[9,3] -2.938 1.089 -5.351 -2.829 -1.128
beta0_pH[10,3] -0.934 0.919 -3.060 -0.806 0.498
beta0_pH[11,3] 0.023 11.675 -22.811 -0.056 23.189
beta0_pH[12,3] -0.087 11.430 -22.179 -0.150 23.085
beta0_pH[13,3] -0.044 11.444 -22.470 -0.043 22.852
beta0_pH[14,3] 0.095 11.594 -22.006 0.143 24.247
beta0_pH[15,3] -0.124 11.337 -22.679 -0.067 21.754
beta0_pH[16,3] -0.238 11.533 -22.876 -0.307 22.403
beta0_pH[11,4] 0.997 0.731 -0.856 1.087 2.166
beta0_pH[12,4] 0.894 0.867 -1.218 0.957 2.383
beta0_pH[13,4] 0.769 0.706 -0.934 0.849 1.966
beta0_pH[14,4] 0.345 0.578 -0.957 0.406 1.305
beta0_pH[15,4] 0.706 0.588 -0.583 0.720 1.840
beta0_pH[16,4] 1.010 0.761 -0.695 1.043 2.384
beta0_pH[11,5] 1.670 0.085 1.523 1.665 1.855
beta0_pH[12,5] -1.409 0.441 -2.377 -1.393 -0.579
beta0_pH[13,5] 2.238 0.255 1.813 2.213 2.809
beta0_pH[14,5] 1.888 0.301 1.386 1.859 2.564
beta0_pH[15,5] 1.111 0.066 0.996 1.106 1.254
beta0_pH[16,5] -0.736 0.055 -0.841 -0.737 -0.631
beta1_pH[1,1] 2.456 0.463 1.755 2.384 3.640
beta1_pH[2,1] 1.592 0.340 1.029 1.567 2.239
beta1_pH[3,1] 1.832 0.491 1.098 1.752 3.167
beta1_pH[4,1] 1.956 0.473 1.238 1.892 3.212
beta1_pH[5,1] 2.334 0.493 1.539 2.286 3.521
beta1_pH[6,1] 2.282 0.870 1.127 2.061 4.433
beta1_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,1] 2.566 0.890 1.285 2.388 4.482
beta1_pH[9,1] 2.173 0.720 1.184 1.994 4.165
beta1_pH[10,1] 2.185 0.392 1.501 2.143 3.015
beta1_pH[11,1] 3.363 0.522 2.619 3.275 4.768
beta1_pH[12,1] 2.888 0.274 2.355 2.883 3.431
beta1_pH[13,1] 3.436 0.250 2.995 3.418 3.979
beta1_pH[14,1] 2.862 0.322 2.268 2.850 3.537
beta1_pH[15,1] 2.444 0.297 1.923 2.419 3.126
beta1_pH[16,1] 4.020 0.311 3.449 4.011 4.674
beta1_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[1,3] 3.458 1.648 1.258 3.019 7.273
beta1_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,3] 4.981 2.399 2.245 4.151 11.439
beta1_pH[9,3] 8.624 2.641 4.406 8.271 14.150
beta1_pH[10,3] 7.315 2.729 3.193 6.767 13.304
beta1_pH[11,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,4] 2.175 0.763 0.781 2.129 3.971
beta1_pH[12,4] 3.280 0.877 1.801 3.200 5.470
beta1_pH[13,4] 1.694 0.748 0.315 1.647 3.427
beta1_pH[14,4] 2.304 0.663 1.126 2.261 3.717
beta1_pH[15,4] 1.746 0.661 0.411 1.744 3.124
beta1_pH[16,4] 3.845 0.816 2.323 3.815 5.655
beta1_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,1] 1.328 0.858 0.341 1.027 3.386
beta2_pH[2,1] 1.607 0.930 0.328 1.385 3.413
beta2_pH[3,1] 1.315 0.847 0.183 1.135 2.934
beta2_pH[4,1] 1.090 0.781 0.204 0.828 2.905
beta2_pH[5,1] 2.147 1.798 0.283 1.671 7.480
beta2_pH[6,1] 0.967 1.586 0.112 0.365 5.971
beta2_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,1] 0.671 1.076 0.124 0.315 4.201
beta2_pH[9,1] 1.110 1.790 0.138 0.477 7.840
beta2_pH[10,1] 1.247 1.344 0.284 0.777 5.449
beta2_pH[11,1] 0.415 0.200 0.181 0.367 0.912
beta2_pH[12,1] 1.581 0.806 0.619 1.341 3.687
beta2_pH[13,1] 1.022 0.437 0.528 0.911 2.198
beta2_pH[14,1] 1.109 0.629 0.427 0.918 3.010
beta2_pH[15,1] 0.995 0.588 0.342 0.846 2.710
beta2_pH[16,1] 1.063 0.331 0.642 0.997 1.870
beta2_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,3] 2.360 4.227 0.113 0.486 15.542
beta2_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,3] 2.443 2.055 0.577 1.515 8.193
beta2_pH[9,3] 0.175 0.092 0.080 0.154 0.382
beta2_pH[10,3] 0.289 0.301 0.097 0.225 0.795
beta2_pH[11,3] -8.434 5.944 -22.569 -7.129 -0.422
beta2_pH[12,3] -8.380 5.906 -22.349 -7.139 -0.446
beta2_pH[13,3] -8.436 5.836 -21.981 -7.246 -0.480
beta2_pH[14,3] -8.475 5.877 -22.323 -7.376 -0.373
beta2_pH[15,3] -8.441 5.851 -22.078 -7.309 -0.494
beta2_pH[16,3] -8.419 5.913 -22.757 -7.364 -0.380
beta2_pH[11,4] -5.406 4.406 -16.679 -3.993 -0.665
beta2_pH[12,4] -5.890 4.355 -17.304 -4.674 -0.994
beta2_pH[13,4] -5.444 4.323 -16.591 -4.051 -0.816
beta2_pH[14,4] -5.582 4.345 -17.003 -4.123 -0.974
beta2_pH[15,4] -5.767 4.372 -16.778 -4.440 -0.874
beta2_pH[16,4] -5.544 4.072 -16.833 -4.193 -1.294
beta2_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta3_pH[1,1] 34.832 1.098 32.764 34.831 37.093
beta3_pH[2,1] 32.856 1.450 30.246 32.771 36.226
beta3_pH[3,1] 33.030 1.589 30.054 32.937 36.745
beta3_pH[4,1] 32.508 1.906 29.134 32.423 36.825
beta3_pH[5,1] 27.570 1.106 25.983 27.458 30.446
beta3_pH[6,1] 33.508 3.588 26.375 33.235 41.653
beta3_pH[7,1] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,1] 35.638 3.745 28.999 35.906 42.331
beta3_pH[9,1] 32.694 3.170 27.905 32.015 40.740
beta3_pH[10,1] 32.629 1.318 30.305 32.562 35.340
beta3_pH[11,1] 31.200 1.133 29.190 31.136 33.483
beta3_pH[12,1] 30.788 0.472 29.842 30.799 31.703
beta3_pH[13,1] 31.438 0.573 30.312 31.425 32.599
beta3_pH[14,1] 32.513 0.721 31.297 32.474 34.003
beta3_pH[15,1] 30.917 0.752 29.380 30.929 32.398
beta3_pH[16,1] 30.131 0.360 29.430 30.129 30.814
beta3_pH[1,2] 34.421 5.493 25.433 34.259 43.493
beta3_pH[2,2] 34.402 5.505 25.472 34.369 43.532
beta3_pH[3,2] 30.370 7.111 18.516 30.451 42.244
beta3_pH[4,2] 30.812 7.262 18.783 31.056 42.345
beta3_pH[5,2] 30.386 7.180 18.654 30.186 42.449
beta3_pH[6,2] 30.467 7.222 18.757 30.435 42.447
beta3_pH[7,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[9,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[10,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[11,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[12,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[13,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[14,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[15,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[16,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[1,3] 39.724 2.671 33.071 40.178 43.692
beta3_pH[2,3] 34.638 5.420 25.560 34.705 43.422
beta3_pH[3,3] 30.579 7.192 18.705 30.608 42.424
beta3_pH[4,3] 30.522 7.285 18.628 30.425 42.437
beta3_pH[5,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[6,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[7,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,3] 42.171 0.826 41.024 41.904 43.878
beta3_pH[9,3] 36.023 1.776 34.524 35.349 41.394
beta3_pH[10,3] 37.805 2.509 34.623 37.195 43.381
beta3_pH[11,3] 38.463 2.597 34.272 38.385 42.785
beta3_pH[12,3] 38.531 2.580 34.250 38.574 42.761
beta3_pH[13,3] 38.466 2.590 34.192 38.483 42.738
beta3_pH[14,3] 38.508 2.584 34.227 38.510 42.776
beta3_pH[15,3] 38.507 2.608 34.222 38.425 42.824
beta3_pH[16,3] 38.487 2.596 34.217 38.536 42.735
beta3_pH[11,4] 40.738 1.213 38.480 40.629 43.229
beta3_pH[12,4] 41.853 0.817 40.240 42.136 42.935
beta3_pH[13,4] 41.429 1.129 39.184 41.391 43.589
beta3_pH[14,4] 40.360 0.743 38.875 40.360 41.765
beta3_pH[15,4] 40.941 0.990 38.556 40.945 42.917
beta3_pH[16,4] 41.735 0.434 40.908 41.787 42.592
beta3_pH[11,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[12,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[13,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[14,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[15,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[16,5] 25.000 0.000 25.000 25.000 25.000
beta4_pH[1,1] 0.263 0.021 0.222 0.263 0.303
beta4_pH[2,1] 0.211 0.005 0.200 0.211 0.221
beta4_pH[3,1] 0.173 0.007 0.159 0.172 0.187
beta4_pH[4,1] 0.293 0.011 0.271 0.294 0.316
beta4_pH[5,1] 0.417 0.262 0.096 0.375 0.938
beta4_pH[6,1] 0.774 0.327 0.319 0.737 1.474
beta4_pH[7,1] 0.866 0.351 0.366 0.812 1.705
beta4_pH[8,1] 0.155 0.068 0.045 0.150 0.301
beta4_pH[9,1] 0.217 0.087 0.062 0.213 0.397
beta4_pH[10,1] 0.001 0.001 0.000 0.001 0.003
beta4_pH[11,1] 0.446 0.012 0.423 0.446 0.468
beta4_pH[12,1] 0.643 0.022 0.600 0.642 0.686
beta4_pH[13,1] 0.631 0.013 0.605 0.631 0.658
beta4_pH[14,1] 0.667 0.030 0.609 0.667 0.729
beta4_pH[15,1] 0.795 0.021 0.755 0.795 0.837
beta4_pH[16,1] 1.215 0.042 1.134 1.215 1.298
beta4_pH[1,2] 0.296 0.097 0.126 0.291 0.504
beta4_pH[2,2] 0.515 0.022 0.473 0.514 0.559
beta4_pH[3,2] 0.446 0.012 0.422 0.446 0.469
beta4_pH[4,2] 0.655 0.037 0.584 0.654 0.729
beta4_pH[5,2] 0.716 0.730 0.104 0.533 2.553
beta4_pH[6,2] 0.922 0.840 0.212 0.750 2.645
beta4_pH[7,2] 0.904 0.477 0.353 0.819 1.939
beta4_pH[8,2] 0.847 0.205 0.486 0.835 1.309
beta4_pH[9,2] 1.144 1.002 0.309 0.879 3.640
beta4_pH[10,2] 0.456 0.051 0.360 0.455 0.562
beta4_pH[11,2] 1.269 0.050 1.172 1.269 1.368
beta4_pH[12,2] 0.608 0.110 0.406 0.603 0.835
beta4_pH[13,2] 0.949 0.044 0.869 0.948 1.038
beta4_pH[14,2] 0.932 0.063 0.812 0.933 1.054
beta4_pH[15,2] 1.009 0.026 0.960 1.010 1.059
beta4_pH[16,2] 1.390 0.074 1.250 1.388 1.540
beta4_pH[1,3] 0.241 0.077 0.104 0.236 0.404
beta4_pH[2,3] 0.050 0.017 0.020 0.048 0.084
beta4_pH[3,3] 0.052 0.010 0.034 0.052 0.071
beta4_pH[4,3] 0.166 0.031 0.108 0.166 0.227
beta4_pH[5,3] 0.760 0.931 0.111 0.553 2.626
beta4_pH[6,3] 0.890 0.765 0.205 0.740 2.525
beta4_pH[7,3] 0.889 0.422 0.336 0.817 1.865
beta4_pH[8,3] 0.389 0.183 0.102 0.365 0.813
beta4_pH[9,3] 0.914 0.766 0.173 0.709 2.713
beta4_pH[10,3] 0.054 0.038 0.003 0.046 0.140
beta4_pH[11,3] 0.615 0.259 0.267 0.570 1.204
beta4_pH[12,3] 0.705 0.228 0.383 0.671 1.242
beta4_pH[13,3] 0.713 0.197 0.406 0.687 1.163
beta4_pH[14,3] 0.860 0.222 0.508 0.834 1.383
beta4_pH[15,3] 0.826 0.252 0.444 0.786 1.427
beta4_pH[16,3] 0.824 0.377 0.320 0.749 1.758
beta4_pH[1,4] 0.636 0.377 0.196 0.553 1.551
beta4_pH[2,4] 0.833 1.022 0.119 0.580 3.035
beta4_pH[3,4] 0.831 0.824 0.113 0.605 2.870
beta4_pH[4,4] 0.734 0.476 0.202 0.621 1.948
beta4_pH[5,4] 0.750 0.766 0.107 0.555 2.556
beta4_pH[6,4] 0.924 0.874 0.214 0.743 2.863
beta4_pH[7,4] 0.892 0.440 0.377 0.808 1.917
beta4_pH[8,4] 0.762 0.678 0.145 0.603 2.303
beta4_pH[9,4] 0.778 0.852 0.127 0.590 2.656
beta4_pH[10,4] 1.457 2.986 0.006 0.368 11.011
beta4_pH[11,4] 0.396 0.016 0.366 0.396 0.427
beta4_pH[12,4] 0.551 0.053 0.450 0.551 0.659
beta4_pH[13,4] 0.834 0.040 0.758 0.834 0.915
beta4_pH[14,4] 0.864 0.072 0.724 0.861 1.006
beta4_pH[15,4] 0.626 0.018 0.593 0.626 0.661
beta4_pH[16,4] 0.651 0.039 0.577 0.649 0.731
beta4_pH[1,5] 0.634 0.369 0.197 0.556 1.543
beta4_pH[2,5] 0.834 0.911 0.102 0.596 3.146
beta4_pH[3,5] 0.860 0.910 0.115 0.598 3.181
beta4_pH[4,5] 0.736 0.497 0.219 0.618 1.984
beta4_pH[5,5] 0.765 0.979 0.107 0.551 2.693
beta4_pH[6,5] 0.925 0.853 0.216 0.731 2.767
beta4_pH[7,5] 0.892 0.418 0.345 0.817 1.919
beta4_pH[8,5] 0.736 0.606 0.149 0.591 2.345
beta4_pH[9,5] 0.777 0.923 0.122 0.561 2.728
beta4_pH[10,5] 1.504 3.134 0.007 0.375 12.287
beta4_pH[11,5] 0.362 0.077 0.215 0.361 0.514
beta4_pH[12,5] 0.693 0.204 0.374 0.663 1.177
beta4_pH[13,5] 0.973 0.112 0.779 0.965 1.216
beta4_pH[14,5] 0.845 0.158 0.570 0.836 1.187
beta4_pH[15,5] 1.304 0.116 1.087 1.301 1.545
beta4_pH[16,5] 0.916 0.376 0.354 0.863 1.843
beta0_pelagic[1] 2.268 0.036 2.202 2.266 2.342
beta0_pelagic[2] 1.537 0.007 1.524 1.537 1.551
beta0_pelagic[3] 0.153 0.285 -0.615 0.214 0.548
beta0_pelagic[4] 0.219 0.230 -0.326 0.238 0.636
beta0_pelagic[5] 1.442 0.157 1.144 1.441 1.749
beta0_pelagic[6] 1.357 0.061 1.239 1.356 1.482
beta0_pelagic[7] 1.470 0.039 1.398 1.469 1.549
beta0_pelagic[8] 1.672 0.044 1.591 1.670 1.763
beta0_pelagic[9] 2.777 0.063 2.662 2.774 2.906
beta0_pelagic[10] 2.604 0.032 2.543 2.603 2.669
beta0_pelagic[11] -0.132 0.154 -0.462 -0.116 0.137
beta0_pelagic[12] 1.800 0.035 1.733 1.798 1.874
beta0_pelagic[13] -1.282 0.678 -2.553 -1.266 -0.140
beta0_pelagic[14] -0.250 0.151 -0.589 -0.238 0.007
beta0_pelagic[15] -0.210 0.035 -0.304 -0.196 -0.179
beta0_pelagic[16] -0.496 0.357 -1.333 -0.420 0.028
beta1_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[3] 0.814 0.314 0.345 0.755 1.655
beta1_pelagic[4] 0.830 0.251 0.388 0.815 1.394
beta1_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[11] 1.150 0.245 0.732 1.132 1.691
beta1_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[13] 3.326 0.919 1.563 3.363 4.857
beta1_pelagic[14] 1.672 0.375 1.132 1.604 2.634
beta1_pelagic[15] 0.096 0.155 0.000 0.016 0.571
beta1_pelagic[16] 1.650 0.548 0.839 1.604 2.870
beta2_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[3] 5.948 5.999 0.168 4.066 20.698
beta2_pelagic[4] 7.412 5.979 0.311 6.269 22.345
beta2_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[11] 1.803 1.506 0.264 1.306 5.945
beta2_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[13] 0.159 0.078 0.075 0.139 0.376
beta2_pelagic[14] 1.016 0.976 0.243 0.699 4.041
beta2_pelagic[15] 2.980 2.226 0.161 2.445 8.352
beta2_pelagic[16] 0.637 0.846 0.121 0.324 3.480
beta3_pelagic[1] 28.644 3.885 22.358 28.450 35.140
beta3_pelagic[2] 28.793 3.923 22.392 28.798 35.154
beta3_pelagic[3] 27.901 2.912 22.927 28.343 33.138
beta3_pelagic[4] 26.102 1.530 22.934 26.274 29.436
beta3_pelagic[5] 29.289 6.659 18.552 29.173 40.337
beta3_pelagic[6] 29.619 6.576 18.638 29.764 40.508
beta3_pelagic[7] 29.422 6.692 18.539 29.385 40.426
beta3_pelagic[8] 29.458 6.672 18.602 29.307 40.387
beta3_pelagic[9] 29.541 6.722 18.436 29.497 40.475
beta3_pelagic[10] 29.523 6.768 18.546 29.433 40.460
beta3_pelagic[11] 30.208 1.125 28.207 30.121 32.723
beta3_pelagic[12] 30.000 0.000 30.000 30.000 30.000
beta3_pelagic[13] 31.044 4.891 23.503 30.544 40.660
beta3_pelagic[14] 33.971 1.490 31.549 33.850 37.501
beta3_pelagic[15] 34.369 5.695 23.717 35.284 41.741
beta3_pelagic[16] 31.290 2.626 26.310 31.176 36.774
beta0_yellow[1] -0.484 0.079 -0.639 -0.483 -0.328
beta0_yellow[2] 0.537 0.087 0.335 0.547 0.674
beta0_yellow[3] -0.218 0.126 -0.466 -0.219 0.010
beta0_yellow[4] 0.872 0.260 0.164 0.935 1.185
beta0_yellow[5] 0.235 0.324 -0.443 0.248 0.813
beta0_yellow[6] 0.819 0.106 0.615 0.817 1.036
beta0_yellow[7] 1.017 0.103 0.823 1.013 1.221
beta0_yellow[8] 1.016 0.094 0.832 1.018 1.195
beta0_yellow[9] 0.442 0.101 0.249 0.443 0.641
beta0_yellow[10] 0.538 0.070 0.405 0.538 0.674
beta0_yellow[11] -0.209 0.130 -0.507 -0.196 -0.002
beta0_yellow[12] -1.087 0.094 -1.267 -1.087 -0.905
beta0_yellow[13] -1.338 0.169 -1.721 -1.325 -1.056
beta0_yellow[14] -0.249 0.153 -0.574 -0.235 -0.003
beta0_yellow[15] -1.245 0.133 -1.564 -1.231 -1.028
beta0_yellow[16] -0.258 0.093 -0.442 -0.259 -0.077
beta1_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[2] 1.021 0.241 0.571 1.011 1.526
beta1_yellow[3] 0.480 0.236 0.026 0.490 0.923
beta1_yellow[4] 1.360 0.630 0.642 1.190 3.259
beta1_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[11] 0.272 0.346 0.001 0.092 1.147
beta1_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[13] 2.112 1.732 0.060 1.324 6.037
beta1_yellow[14] 0.694 0.916 0.013 0.433 3.513
beta1_yellow[15] 0.163 0.214 0.011 0.075 0.737
beta1_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[2] -7.655 5.874 -22.044 -6.391 -0.367
beta2_yellow[3] -7.707 5.872 -21.543 -6.468 -0.491
beta2_yellow[4] -5.353 5.778 -19.777 -3.195 -0.107
beta2_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[11] -7.690 5.820 -22.276 -6.528 -0.308
beta2_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[13] -4.831 4.687 -16.709 -3.201 -0.318
beta2_yellow[14] -6.235 5.092 -19.024 -4.935 -0.354
beta2_yellow[15] -7.665 5.592 -21.038 -6.500 -0.333
beta2_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta3_yellow[1] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[2] 28.927 1.195 27.300 28.678 32.111
beta3_yellow[3] 32.513 3.068 24.106 32.630 38.415
beta3_yellow[4] 30.050 3.315 23.932 29.061 35.670
beta3_yellow[5] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[6] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[7] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[8] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[9] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[10] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[11] 35.377 3.866 30.378 34.294 43.546
beta3_yellow[12] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[13] 32.096 2.679 29.079 31.584 37.793
beta3_yellow[14] 33.855 3.536 29.211 33.392 42.854
beta3_yellow[15] 36.576 3.796 30.368 36.408 43.577
beta3_yellow[16] 35.000 0.000 35.000 35.000 35.000
beta0_black[1] 0.014 0.077 -0.134 0.015 0.167
beta0_black[2] 1.855 0.032 1.793 1.854 1.918
beta0_black[3] 1.240 0.037 1.168 1.239 1.313
beta0_black[4] 2.426 0.047 2.337 2.426 2.520
beta0_black[5] 5.132 2.361 1.887 4.685 11.030
beta0_black[6] 5.267 2.214 2.607 4.659 11.180
beta0_black[7] 4.119 2.049 1.639 3.601 9.689
beta0_black[8] 0.885 0.217 0.476 0.874 1.356
beta0_black[9] 2.813 0.272 2.301 2.804 3.366
beta0_black[10] 1.287 0.051 1.189 1.287 1.385
beta0_black[11] 2.978 0.051 2.879 2.978 3.079
beta0_black[12] 4.920 0.191 4.549 4.916 5.289
beta0_black[13] 0.417 0.322 0.017 0.340 1.220
beta0_black[14] 2.850 0.111 2.628 2.850 3.070
beta0_black[15] 0.722 0.051 0.619 0.723 0.820
beta0_black[16] 3.773 0.073 3.629 3.773 3.913
beta2_black[1] 0.000 0.000 0.000 0.000 0.000
beta2_black[2] 0.000 0.000 0.000 0.000 0.000
beta2_black[3] 0.000 0.000 0.000 0.000 0.000
beta2_black[4] 0.000 0.000 0.000 0.000 0.000
beta2_black[5] 0.000 0.000 0.000 0.000 0.000
beta2_black[6] 0.000 0.000 0.000 0.000 0.000
beta2_black[7] 0.000 0.000 0.000 0.000 0.000
beta2_black[8] 0.000 0.000 0.000 0.000 0.000
beta2_black[9] 0.000 0.000 0.000 0.000 0.000
beta2_black[10] 0.000 0.000 0.000 0.000 0.000
beta2_black[11] 0.000 0.000 0.000 0.000 0.000
beta2_black[12] 0.000 0.000 0.000 0.000 0.000
beta2_black[13] -2.590 1.769 -6.882 -2.242 -0.325
beta2_black[14] 0.000 0.000 0.000 0.000 0.000
beta2_black[15] 0.000 0.000 0.000 0.000 0.000
beta2_black[16] 0.000 0.000 0.000 0.000 0.000
beta3_black[1] 25.000 0.000 25.000 25.000 25.000
beta3_black[2] 25.000 0.000 25.000 25.000 25.000
beta3_black[3] 25.000 0.000 25.000 25.000 25.000
beta3_black[4] 25.000 0.000 25.000 25.000 25.000
beta3_black[5] 25.000 0.000 25.000 25.000 25.000
beta3_black[6] 25.000 0.000 25.000 25.000 25.000
beta3_black[7] 25.000 0.000 25.000 25.000 25.000
beta3_black[8] 25.000 0.000 25.000 25.000 25.000
beta3_black[9] 25.000 0.000 25.000 25.000 25.000
beta3_black[10] 25.000 0.000 25.000 25.000 25.000
beta3_black[11] 25.000 0.000 25.000 25.000 25.000
beta3_black[12] 25.000 0.000 25.000 25.000 25.000
beta3_black[13] 37.453 3.339 26.008 38.184 40.950
beta3_black[14] 25.000 0.000 25.000 25.000 25.000
beta3_black[15] 25.000 0.000 25.000 25.000 25.000
beta3_black[16] 25.000 0.000 25.000 25.000 25.000
beta4_black[1] -0.329 0.092 -0.509 -0.329 -0.145
beta4_black[2] 0.260 0.053 0.156 0.260 0.364
beta4_black[3] -0.722 0.064 -0.844 -0.721 -0.595
beta4_black[4] 0.630 0.131 0.386 0.629 0.889
beta4_black[5] 0.770 1.483 -1.242 0.499 4.445
beta4_black[6] 0.887 1.606 -0.794 0.559 4.230
beta4_black[7] 0.768 1.412 -0.921 0.485 4.315
beta4_black[8] 0.080 0.290 -0.515 0.083 0.627
beta4_black[9] 1.022 0.797 -0.121 0.883 2.994
beta4_black[10] 0.073 0.067 -0.053 0.073 0.205
beta4_black[11] -0.607 0.065 -0.734 -0.605 -0.484
beta4_black[12] 1.232 0.290 0.694 1.219 1.834
beta4_black[13] -1.004 0.077 -1.163 -1.001 -0.862
beta4_black[14] 0.229 0.180 -0.112 0.227 0.593
beta4_black[15] 0.003 0.046 -0.085 0.002 0.095
beta4_black[16] -0.416 0.117 -0.639 -0.421 -0.181
beta0_dsr[11] -0.337 0.009 -0.356 -0.337 -0.319
beta0_dsr[12] 1.421 0.044 1.336 1.421 1.508
beta0_dsr[13] 0.130 0.021 0.088 0.131 0.173
beta0_dsr[14] -0.135 0.018 -0.168 -0.135 -0.100
beta0_dsr[15] 0.184 0.014 0.157 0.184 0.212
beta0_dsr[16] 0.208 0.015 0.178 0.208 0.238
beta1_dsr[11] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[12] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[13] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[14] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[15] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[16] 0.000 0.000 0.000 0.000 0.000
beta4_dsr[11] 0.549 0.035 0.482 0.548 0.620
beta4_dsr[12] 0.818 0.239 0.354 0.811 1.285
beta4_dsr[13] -0.046 0.041 -0.126 -0.047 0.032
beta4_dsr[14] -0.206 0.071 -0.345 -0.206 -0.068
beta4_dsr[15] 0.237 0.022 0.194 0.237 0.280
beta4_dsr[16] -0.100 0.046 -0.190 -0.100 -0.011
beta0_slope[11] -2.200 0.021 -2.243 -2.201 -2.158
beta0_slope[12] -5.429 0.225 -5.873 -5.428 -4.985
beta0_slope[13] -1.456 0.034 -1.526 -1.456 -1.390
beta0_slope[14] -2.693 0.051 -2.797 -2.691 -2.599
beta0_slope[15] -1.566 0.022 -1.610 -1.567 -1.524
beta0_slope[16] -3.013 0.047 -3.106 -3.012 -2.924
beta1_slope[11] 0.000 0.000 0.000 0.000 0.000
beta1_slope[12] 0.000 0.000 0.000 0.000 0.000
beta1_slope[13] 0.000 0.000 0.000 0.000 0.000
beta1_slope[14] 0.000 0.000 0.000 0.000 0.000
beta1_slope[15] 0.000 0.000 0.000 0.000 0.000
beta1_slope[16] 0.000 0.000 0.000 0.000 0.000
beta4_slope[11] -0.342 0.077 -0.497 -0.339 -0.194
beta4_slope[12] -4.673 0.812 -6.413 -4.621 -3.182
beta4_slope[13] 0.281 0.054 0.172 0.282 0.384
beta4_slope[14] -0.182 0.155 -0.495 -0.178 0.120
beta4_slope[15] -0.387 0.034 -0.454 -0.387 -0.319
beta4_slope[16] 0.218 0.104 0.008 0.219 0.410
beta5_slope[11] 3.081 0.108 2.869 3.077 3.294
beta5_slope[12] 2.006 0.389 1.278 1.995 2.788
beta5_slope[13] 0.934 0.081 0.776 0.934 1.093
beta5_slope[14] 3.800 0.362 3.140 3.778 4.565
beta5_slope[15] 1.388 0.071 1.253 1.388 1.528
beta5_slope[16] 2.679 0.234 2.245 2.671 3.172
tau_comp[1] 3.784 0.294 3.245 3.773 4.392
tau_comp[2] 2.854 0.257 2.381 2.842 3.394
tau_comp[3] 0.867 0.073 0.730 0.864 1.018
tau_comp[4] 1.279 0.195 0.941 1.268 1.698
tau_comp[5] 1.201 0.180 0.889 1.190 1.589
sigma_H[1] 0.205 0.049 0.118 0.202 0.311
sigma_H[2] 0.175 0.028 0.127 0.173 0.236
sigma_H[3] 0.187 0.041 0.114 0.184 0.273
sigma_H[4] 0.464 0.081 0.334 0.457 0.645
sigma_H[5] 0.964 0.202 0.604 0.954 1.394
sigma_H[6] 0.535 0.192 0.162 0.527 0.923
sigma_H[7] 0.310 0.062 0.212 0.300 0.456
sigma_H[8] 0.416 0.089 0.265 0.408 0.604
sigma_H[9] 0.512 0.117 0.335 0.498 0.777
sigma_H[10] 0.208 0.039 0.142 0.205 0.290
sigma_H[11] 0.267 0.044 0.194 0.262 0.366
sigma_H[12] 0.388 0.152 0.200 0.342 0.740
sigma_H[13] 0.218 0.036 0.157 0.215 0.298
sigma_H[14] 0.499 0.090 0.340 0.493 0.685
sigma_H[15] 0.235 0.038 0.170 0.233 0.316
sigma_H[16] 0.211 0.040 0.144 0.206 0.301
lambda_H[1] 2.173 3.667 0.082 0.994 11.315
lambda_H[2] 5.309 4.967 0.439 3.816 19.506
lambda_H[3] 5.173 7.245 0.254 2.795 23.527
lambda_H[4] 0.006 0.004 0.001 0.005 0.018
lambda_H[5] 2.835 7.293 0.027 0.692 20.859
lambda_H[6] 6.338 14.399 0.008 0.296 47.364
lambda_H[7] 0.012 0.009 0.002 0.010 0.034
lambda_H[8] 7.908 10.000 0.102 4.319 37.939
lambda_H[9] 0.015 0.010 0.003 0.013 0.042
lambda_H[10] 0.222 0.340 0.027 0.149 0.784
lambda_H[11] 0.235 0.377 0.011 0.109 1.144
lambda_H[12] 3.664 4.938 0.175 2.134 16.835
lambda_H[13] 3.077 2.782 0.274 2.321 10.513
lambda_H[14] 3.334 3.975 0.200 2.155 13.788
lambda_H[15] 0.019 0.022 0.003 0.014 0.061
lambda_H[16] 0.507 0.805 0.035 0.260 2.628
mu_lambda_H[1] 4.090 1.890 1.061 3.893 8.232
mu_lambda_H[2] 3.698 1.954 0.518 3.516 7.816
mu_lambda_H[3] 3.354 1.834 0.690 3.062 7.521
sigma_lambda_H[1] 8.184 4.347 1.684 7.513 17.831
sigma_lambda_H[2] 8.132 4.699 0.851 7.543 18.165
sigma_lambda_H[3] 6.199 4.067 0.941 5.354 16.294
beta_H[1,1] 6.459 1.345 3.206 6.673 8.437
beta_H[2,1] 9.763 0.602 8.394 9.802 10.824
beta_H[3,1] 7.827 0.801 5.940 7.928 9.125
beta_H[4,1] 10.629 7.814 -4.911 10.668 25.512
beta_H[5,1] -0.326 2.491 -5.531 -0.194 3.956
beta_H[6,1] 3.184 4.310 -7.329 4.596 8.379
beta_H[7,1] 1.335 5.648 -10.704 1.721 11.530
beta_H[8,1] 1.333 3.522 -2.380 1.251 3.571
beta_H[9,1] 13.330 5.731 2.377 13.142 25.120
beta_H[10,1] 6.837 1.813 3.005 6.896 10.338
beta_H[11,1] 4.888 3.533 -3.075 5.546 9.903
beta_H[12,1] 2.583 1.109 0.690 2.507 4.979
beta_H[13,1] 9.091 0.884 7.262 9.164 10.552
beta_H[14,1] 2.129 1.033 0.155 2.120 4.267
beta_H[15,1] -7.014 3.468 -13.613 -7.103 0.402
beta_H[16,1] 3.934 3.025 -0.844 3.550 10.785
beta_H[1,2] 7.781 0.273 7.201 7.786 8.295
beta_H[2,2] 9.884 0.167 9.560 9.881 10.230
beta_H[3,2] 8.799 0.202 8.406 8.795 9.205
beta_H[4,2] 3.007 1.468 0.321 2.988 6.195
beta_H[5,2] 1.775 0.994 -0.178 1.787 3.729
beta_H[6,2] 6.052 1.106 3.415 6.225 7.801
beta_H[7,2] 2.376 1.059 0.396 2.309 4.652
beta_H[8,2] 3.009 1.083 1.254 3.150 4.257
beta_H[9,2] 3.222 1.102 0.942 3.217 5.459
beta_H[10,2] 8.035 0.362 7.294 8.047 8.706
beta_H[11,2] 9.820 0.639 8.860 9.725 11.252
beta_H[12,2] 3.934 0.380 3.227 3.922 4.747
beta_H[13,2] 9.152 0.243 8.692 9.142 9.662
beta_H[14,2] 4.030 0.344 3.352 4.023 4.708
beta_H[15,2] 11.532 0.623 10.219 11.543 12.738
beta_H[16,2] 4.467 0.818 2.990 4.437 6.133
beta_H[1,3] 8.456 0.263 7.976 8.442 8.993
beta_H[2,3] 9.855 0.153 9.555 9.853 10.164
beta_H[3,3] 9.464 0.160 9.171 9.458 9.793
beta_H[4,3] -2.185 0.876 -3.959 -2.178 -0.462
beta_H[5,3] 3.895 0.612 2.681 3.909 5.090
beta_H[6,3] 8.408 1.185 6.776 8.088 10.884
beta_H[7,3] -2.557 0.633 -3.858 -2.550 -1.342
beta_H[8,3] 5.300 0.513 4.684 5.228 6.449
beta_H[9,3] -2.733 0.707 -4.158 -2.729 -1.358
beta_H[10,3] 8.545 0.275 8.017 8.537 9.113
beta_H[11,3] 8.494 0.292 7.863 8.513 8.994
beta_H[12,3] 5.264 0.325 4.513 5.311 5.792
beta_H[13,3] 8.862 0.175 8.504 8.867 9.195
beta_H[14,3] 5.785 0.268 5.196 5.802 6.265
beta_H[15,3] 10.282 0.295 9.708 10.286 10.859
beta_H[16,3] 5.946 0.594 4.704 5.976 7.012
beta_H[1,4] 8.093 0.186 7.689 8.107 8.416
beta_H[2,4] 9.852 0.127 9.588 9.856 10.089
beta_H[3,4] 9.954 0.153 9.613 9.964 10.222
beta_H[4,4] 11.808 0.468 10.886 11.810 12.723
beta_H[5,4] 5.559 0.735 4.374 5.480 7.182
beta_H[6,4] 7.219 0.923 5.297 7.452 8.491
beta_H[7,4] 8.545 0.347 7.878 8.542 9.223
beta_H[8,4] 6.726 0.261 6.182 6.742 7.138
beta_H[9,4] 7.426 0.447 6.564 7.419 8.316
beta_H[10,4] 7.519 0.219 7.099 7.513 7.960
beta_H[11,4] 9.504 0.194 9.119 9.505 9.887
beta_H[12,4] 7.249 0.205 6.842 7.247 7.679
beta_H[13,4] 9.062 0.139 8.787 9.063 9.327
beta_H[14,4] 7.832 0.215 7.420 7.826 8.283
beta_H[15,4] 9.462 0.220 9.031 9.462 9.881
beta_H[16,4] 9.434 0.236 9.013 9.420 9.912
beta_H[1,5] 8.900 0.153 8.596 8.902 9.187
beta_H[2,5] 10.673 0.098 10.491 10.669 10.878
beta_H[3,5] 10.764 0.158 10.479 10.759 11.090
beta_H[4,5] 8.356 0.497 7.394 8.351 9.356
beta_H[5,5] 5.201 0.595 3.786 5.279 6.155
beta_H[6,5] 9.010 0.640 8.055 8.872 10.438
beta_H[7,5] 6.575 0.346 5.894 6.574 7.259
beta_H[8,5] 8.253 0.218 7.906 8.233 8.700
beta_H[9,5] 8.190 0.458 7.306 8.195 9.078
beta_H[10,5] 10.042 0.228 9.608 10.044 10.503
beta_H[11,5] 11.592 0.225 11.165 11.590 12.038
beta_H[12,5] 8.613 0.192 8.232 8.613 8.991
beta_H[13,5] 10.135 0.129 9.886 10.132 10.389
beta_H[14,5] 9.246 0.216 8.840 9.236 9.693
beta_H[15,5] 11.314 0.236 10.866 11.311 11.784
beta_H[16,5] 9.861 0.181 9.487 9.870 10.192
beta_H[1,6] 10.124 0.240 9.748 10.091 10.717
beta_H[2,6] 11.471 0.116 11.244 11.472 11.702
beta_H[3,6] 10.616 0.155 10.287 10.624 10.901
beta_H[4,6] 12.892 0.893 11.075 12.904 14.625
beta_H[5,6] 5.681 0.638 4.532 5.650 6.957
beta_H[6,6] 8.741 0.589 7.401 8.777 9.802
beta_H[7,6] 10.548 0.588 9.370 10.549 11.708
beta_H[8,6] 9.530 0.281 8.992 9.552 9.961
beta_H[9,6] 8.298 0.755 6.810 8.288 9.818
beta_H[10,6] 9.140 0.349 8.399 9.159 9.764
beta_H[11,6] 10.684 0.366 9.911 10.703 11.337
beta_H[12,6] 9.377 0.240 8.920 9.369 9.891
beta_H[13,6] 11.092 0.152 10.800 11.092 11.399
beta_H[14,6] 9.923 0.281 9.344 9.925 10.472
beta_H[15,6] 10.515 0.424 9.661 10.529 11.342
beta_H[16,6] 10.587 0.254 10.041 10.594 11.067
beta_H[1,7] 10.153 1.272 6.848 10.436 11.933
beta_H[2,7] 12.116 0.515 10.971 12.146 13.061
beta_H[3,7] 10.286 0.676 8.841 10.343 11.492
beta_H[4,7] 1.530 4.648 -7.715 1.545 10.630
beta_H[5,7] 5.819 2.123 1.926 5.906 9.819
beta_H[6,7] 7.180 3.226 -1.162 8.285 11.198
beta_H[7,7] 7.644 2.933 1.923 7.612 13.436
beta_H[8,7] 10.913 1.017 9.344 10.881 12.578
beta_H[9,7] 4.862 3.934 -2.806 4.819 12.789
beta_H[10,7] 10.035 1.654 7.094 9.883 13.548
beta_H[11,7] 11.397 1.824 8.087 11.215 15.370
beta_H[12,7] 9.782 0.978 7.654 9.836 11.464
beta_H[13,7] 11.705 0.699 10.101 11.780 12.871
beta_H[14,7] 10.624 0.944 8.784 10.654 12.382
beta_H[15,7] 14.099 2.213 9.872 14.093 18.547
beta_H[16,7] 12.886 1.405 10.427 12.752 15.917
beta0_H[1] 8.758 18.205 -26.678 8.527 47.247
beta0_H[2] 10.477 7.948 -6.137 10.643 25.570
beta0_H[3] 9.718 10.420 -9.657 9.561 31.094
beta0_H[4] 9.281 183.032 -352.262 4.684 376.742
beta0_H[5] 4.866 27.150 -47.308 4.697 58.321
beta0_H[6] 7.011 51.697 -106.590 7.988 116.911
beta0_H[7] 4.310 138.373 -288.386 7.318 277.377
beta0_H[8] 5.922 34.588 -18.495 6.467 29.430
beta0_H[9] 8.567 121.396 -226.680 8.217 251.448
beta0_H[10] 7.458 37.724 -71.313 8.188 78.673
beta0_H[11] 9.537 50.242 -95.678 9.071 115.746
beta0_H[12] 6.529 12.340 -17.362 6.594 30.900
beta0_H[13] 9.564 10.621 -12.123 9.863 28.797
beta0_H[14] 7.041 12.252 -17.503 7.317 30.363
beta0_H[15] 8.653 111.546 -208.780 8.225 238.551
beta0_H[16] 7.390 30.665 -55.148 6.992 73.691